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Abstract 

Siegert pseudostates are purely outgoing states at some fixed point expanded 
over a finite basis. With discretized variables, they provide an accurate de- 
scription of scattering in the s wave for short-range potentials with few basis 
states. The i?-matrix method combined with a Lagrange basis, i.e. functions 
which vanish at all points of a mesh but one, leads to simple mesh-like equa- 
tions which also allow an accurate description of scattering. These methods 
are shown to be exactly equivalent for any basis size, with or without dis- 
cretization. The comparison of their assumptions shows how to accurately 
derive poles of the scattering matrix in the i?-matrix formalism and suggests 
how to extend the Siegert-pseudostate method to higher partial waves. The 
different concepts are illustrated with the Bargmann potential and with the 
centrifugal potential. A simplification of the ii-matrix treatment can usefully 
be extended to the Siegert-pseudostate method. 

PACS numbers: 03.65.Nk, ll.55.-m, 02.70.Hm 
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I. INTRODUCTION 



Recently a new approach for the description of scattering using Siegert pseudostates has 
been proposed |],|2| . Siegert states are bounded solutions of the Schrodinger equation which 
are purely outgoing at infinity JJ]. These states are particularly interesting because their 
complex wave numbers provide the poles of the scattering matrix. However the derivation of 
these states and their use is not computationally simple. Therefore, Tolstikhin, Ostrovsky, 
and Nakamura [|I],|2| have proposed to use modified Siegert states which satisfy purely out- 
going conditions at some finite distance. These modified states correspond to exact Siegert 
states of a truncated potential. The cutoff distance is similar to the channel radius in the 
.R-matrix theory and its introduction opens a way to possible relations between the two 
methods. 

The Siegert pseudostates are defined by an expansion of the Siegert states of the trun- 
cated potential over a finite basis which becomes complete when the number of basis states 
tends to infinity. In Refs. the authors present an efficient way for deriving the Siegert 
pseudostates in the s wave for short-range potentials. They establish a number of remark- 
able mathematical properties of these pseudostates and of the corresponding complex wave 
numbers. This method can be simplified with the help of a Gauss quadrature in the spirit 
of the discrete- variable representation Then a simple matrix representation is obtained 
where the potential matrix is diagonal. 

The i?-matrix theory is a powerful tool, not only to parametrize scattering matrices and 
cross sections but also to solve the Schrodinger equation at positive energies 0-0. In this 
method, the configuration space is divided in two parts, separated at the channel radius. In 
the external part, the wave functions are approximated by their asymptotic expressions. In 
the internal part, a finite basis of square- integrable functions can be used. 

The Lagrange-mesh method is an approximate variational calculation which resembles a 



mesh calculation pHT0|. This property is obtained by using a basis of Lagrange functions, 
i.e. orthonormal functions which vanish at all points of an associated mesh but one, and the 
Gauss quadrature corresponding to this mesh. In spite of its simplicity, the accuracy of the 
Lagrange-mesh method is very high, a property not explained yet [|TO|] . 

In the single-channel case, Malegat JT1J combined the i?-matrix theory with a Lagrange- 
mesh method based on shifted Legendre polynomials to study the scattering by a simple 
solvable potential. Strikingly, the accuracy of the .R-matrix method on a mesh is as good 
as the accuracy of the i?-matrix method using the corresponding Lagrange basis without 
any approximation [12j . This method can easily be extended to multichannel scattering and 
gives accurate results for realistic problems [|13[ . 



In the present paper, we show that the Siegert-pseudostate method and the .R-matrix 
method on a Lagrange mesh are completely equivalent for any basis size when the bases used 
in both approaches are identical. More strikingly both methods remain exactly equivalent 
when their respective mesh approximations are employed. This equivalence sheds new lights 
on both approaches. It emphasizes the poorly known fact that the .R-matrix method can give 
a direct access to the poles of the scattering matrix. The technique presented in Refs. 
provides an accurate practical way of solving this problem for the s wave of short-range 
potentials. Symetrically, the validity of the .R-matrix method for higher partial waves and for 
long-range potentials indicates the way to natural generalizations of the Siegert-pseudostate 
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method. 

The i?-matrix method is summarized in Sec. II. Its application to the Lagrange- Legendre 
mesh is presented in Sec. III. The equivalence with the Siegert-pseudostate method is dis- 
cussed in Sec. IV. Examples are commented in Sec. V. Concluding remarks are presented in 
Sec. VI. 

II. i?-MATRIX METHOD 

As in Ref. [[|, we restrict ourselves to a single channel. Contrary to that reference, we 
first consider an arbitrary partial wave and potentials with a possible Coulomb asymptotic 
behavior. We follow the notations of Ref. ||12|| . A translation into the notations of Ref. |2| 
is delayed to Sec. IV (see also Table | in that section). 

For the Zth partial wave, the radial Schrodinger equation can be written as 

(H t - E)u h = (1) 
with ui(0) = 0. With fi — m — 1, the radial Hamiltonian reads 

* = * + ™ = ^-£ + ^) + ™, < 2 > 

where V(r) is a radial potential. The phase shift Si is obtained from the asymptotic behavior 
of bounded solutions 



u i( r ) — > k{kr) - SiOi(kr), (3) 

r— >oo 

for positive values of the wave number k corresponding to the energy E = \k 2 . The functions 
Ii and 0\ are ingoing and outgoing Coulomb functions and Si = exp(2i5i) is the scattering 
matrix. In the following k will also take complex values. 

In the i?-matrix method, the configuration space is divided at the channel radius a into 
an internal region and an external region. In the external region, the wave function is 
approximated by its asymptotic form (^). In the internal region, it is expanded on some 
basis. The formalism is conveniently expressed with the help of the Bloch surface operator 
U 

£(fl)= l f(r _ „)(-_!), (4) 

where B is the boundary parameter. The Bloch-Schrodinger equation reads 

(Ht + C(B) - E)ui = C(B)u h (5) 

where the operator Hi + C(B) is Hermitian when B is real. The approximation consists 
in using the asymptotic form (|3]) in the right-hand side of Eq. ([|). The main advantage of 
the i?-matrix method is that an expansion in square-integrable functions can be used in the 
left-hand side. 
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Let us consider a set of N basis functions fi(r) (not necessarily orthogonal) and let us 
expand ui in the internal region as 

N 

M r ) = £c J -/ i (r). ( 6 ) 

3=1 

Equation (||) becomes after projection on fi(r), 

N 

E [Ca{B) - EN l3 }c 3 = (MCiBM - SA>. (7) 

3 = 1 

The elements of the symmetric matrices C and N are defined as 

C lJ (B) = (f l \T l + C(B) + V\f J ) 

= Cy(0) - (BMf^f^a) (8) 

and 

% = (h\f 3 ). (9) 

They correspond to one- dimensional integrals over the variable r from to a. Matrix N 
reduces to the unit matrix when the basis is orthonormal. The right-hand side of Eq. (0) is 
even simpler. Because of the Bloch operator, it only involves values at r = a. 

The coefficients Cj are obtained by solving Eq. (0). The continuity condition at r = a 
between the internal approximation (||) and the asymptotic expression (§) reads 

TV 

Y.c 3 fM) = h(ka) - SMka). (10) 

3=1 

Let us introduce the external logarithmic derivative Li at the channel radius 
and the dimensionless R matrix 

N 

MB) = (2a)- 1 J2 fi{a)[C{B) - EN^f^a). (12) 

i,3=l 

The dependence of the R matrix on the energy E is implied. Introducing the Cj in Eq. (|10"D 
and using Eqs. ( |T2"D and (|TTD , one obtains the scattering matrix for the Ith partial wave 

_ LiM) l-m-B)R l {B) 
1 Oi{ka) l-{L l -B)R l {BY 1 ' 

where L\ is the conjugate of L\. Expression ([13]) has the striking property that it does not 
depend on the boundary parameter B, independently of the size of the basis. Indeed, from 
the matrix relation in Appendix A, one obtains 
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m=RM +B (14) 

for any B, real or complex. Introducing relation ( |I4] ) in Eq. (p~3| ) shows that any B value 
leads to the same scattering matrix as for B = 0. Equation (|14]) is well known in i?-matrix 
theory (see Eq. (IV.2.5a) of Ref. 0). However its validity for the approximation (|T2|) for 
any basis size |l5j is sometimes overlooked. 

The wave function in the internal region is then given by 

ui(r) = fiaR^B^ihika) - S^ka)) 

N N 

xE/.WEM-^/.W. (is) 

3=1 i=l 

As the scattering matrix Si and the external wave function Ii(kr) — SiOi(kr), this expression 
does not depend on the choice for B. Indeed, with the help of relation (|X3|), one easily 
shows that, for any B, it is equal to the similar expression where B is replaced by zero. 



In Refs. [12,13 , the parameter B was chosen equal to zero for obvious reasons of sim- 



plicity. Another interesting choice is [|!4l ,|l5| 



B = L h (16) 

This complex value leads to a complex function Ri(Li) which is not an R matrix in the strict 
sense since R matrices are real. However it is also given by expression (|T^). Equation ( |T3"D 



then takes the simpler form [15 



Si = e" 2 ^^ [1 + 2iP l (ka)R l (L l )l (17) 
where Pi(ka) is defined as 

Pi = i(U-L!) (18) 

and (pi(ka) is half the phase of Oi(ka)/Ii(ka). When k is real, Pi is the penetration factor 
given by the imaginary part of Li and <pi(ka) is the hard-sphere phase shift ||. Since Eq. (|17D 
has no denominator, a direct relation appears between the poles of the scattering matrix 
and of the complex R matrix. 

III. R MATRIX ON A LAGRANGE MESH 

The previous section is valid for arbitrary bases. The calculation of the elements of 
matrix C, 

C l3 {B) = (m + £(0)|/,) - (S/2a)/ i (a)/,(a) + (fi\V\f)), (19) 

involves an evaluation of the matrix elements of the potential, which can be tedious and 
must be repeated when the potential changes. By chosing a Lagrange basis and using the 



associated Gauss quadrature, one can avoid this calculation without losing accuracy [12 



5 



As Lagrange basis, we use functions based on Legendre polynomials |]8|JTT|-[T3| . These 
functions are denoted as /, in Ref. [12]. Here as in Ref. [13] we drop the 'hat' because we 
shall not use any other basis. A Lagrange basis is a set of N functions fi(x) associated with 
a Lagrange mesh of N points axi on the interval [0, a] ||||. The x^s are zeros of the shifted 
Legendre polynomial P/v(2x — 1) [O, i.e., 



P N (2 Xi - 1) = 0. (20) 
The Lagrange functions are continuous and indefinitely differentiable anywhere. They read 



Mr) = (-ir-a 



a-/ .-1/2 rP N [2(r/a) - 1] 



x , 



r — axj 



They satisfy the Lagrange conditions 

fi(axj) = (aXiy 1/2 5ij, 



(21) 



(22) 



i.e., each fi vanishes at all mesh points axj, except at axi. The coefficients Aj are the 
weights associated with a Gauss- Legendre quadrature approximation for the [0, 1] interval. 
The Gauss quadrature on the [0, a] interval reads fl~6 



N 



g(r)dr ps a^\ k g(ax k ). 



(23) 



k=l 



The weights Aj are equal to the traditional Gauss-Legendre weights for the [—1, +1] interval, 
divided by 2. 

The Lagrange functions (|2l|) are not orthogonal [[12] 



(/i|/i> = Sij + 



1 



2N + 1\ 



(l-Xi^l-Xj) 



X 2 j 



(24) 



Because of the Lagrange conditions fl22|) , they are approximately orthogonal at the Gauss 
approximation (E3), 



/, I , i Gauss j- 

\Ji\Jj) = d ir 



(25) 



Strikingly, using the Gauss approximation does not seem to reduce the accuracy of the 
i?-matrix method 1112 . 



At the Gauss approximation, the potential matrix 

(Wlf^^Viax^ (26) 

is diagonal because of Eq. (|22|), and easy to compute. The other matrix elements are exactly 
calculated with the Gauss quadrature. The matrix elements of the sum of the radial part of 
the kinetic energy and of the Bloch operator are given by 



6a 2 Xi(l — Xi) 



AN(N + 1) + 3 + 



1 — Qxi 

X-i ( 1 Xj) 



(27) 
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and, for i ^ j, 



(fi\T +£(p)\fj) 



2a 2 [x i x j {l - Xi){l - xj)} 1 / 2 



N(N + !) + ! + 



(28) 



Finally, the remaining necessary expressions read 

(f t \r- 2 \f,) = a~ 2 x- 2 6 lJ (29) 

and 

(U\r- l 5{r - a)\fj) = a' 2 (-1)^1x^(1 - x<)(l - x,)]^ 2 . (30) 



IV. EQUIVALENCE OF THE SIEGERT-PSEUDOSTATE AND i?-MATRIX 

METHODS 

A. The Siegert-pseudostate method 

First we briefly summarize the Siegert-pseudostate method. For this, we keep the R- 
matrix notations introduced above. Also, we start with definitions for an arbitrary partial 
wave. To avoid confusions, the equations of Ref. [[| are denoted as the reference number 
followed by the equation number in that reference. 

The Siegert pseudostates (j>i n \r) with complex wave numbers ki are solutions of the 
equation 

(^ + £(LW)-^ (n) )# ) = 0, (31) 

where is calculated for k = fcj and = ^k^ 2 . Indeed, with the choice B = Li, the 
right-hand side of Eq. (|7|) vanishes for purely outgoing waves such as the Siegert pseudostates. 
For I = and short-range potentials, definition fl3"I|) is exactly equivalent to the pair of 
equations ||-(la) and @-(lc') (see Ref. [0)- Thanks to the use of the Bloch operator, 
Eq. ( |3~TD is more compact. 



Expanding (f)^ in the internal region as 



N 

*« W (r) = £c< B) /i(r), (32) 

3=1 

one obtains the homogeneous part of Eq. (0), 

TV 

J2lCv(Li n) )-El n) N tJ ]c^ = 0. (33) 
i=i 

However, the dependence of this equation on its eigenvalues fcj is strongly non linear 
because k appears not only in E = \k 2 but also in L\. 
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We now specialize to I = and short-range potentials. Then the free outgoing wave is 
given by Oq = exp(ifcr) and the logarithmic derivative reads 



With the help of Eq. 



L (ka) = ika. 
the system (|33"D can be written as 



N 



£ [Ck(0) - |<C/i(a)/i( 



l fc (n)2 



2^0 



(«) 



0. 



(34) 



(35) 



The eigenvalues appear linearly and quadratically because of the simple form ( |3~4"D of 
L . Ref. || provides an efficient algebraic algorithm for solving the system (|35D , which 
leads to 2iV eigenstates. The quadratic matrix eigenvalue problem is replaced by a standard 
generalized eigenvalue problem of double size. Such a simple algorithm is not yet available 
for system (|33|) with / > 0. 

The fact that Eq. (|3~5D can be solved algebraically gives an access to the physical poles 
of the scattering matrix related to the true Siegert states. The other obtained poles are 
either Siegert states of the truncated potential or non-converged Siegert states (see example 
below). Tolstikhin, Ostrovsky, and Nakamura have shown [Eq. [@]-(59)] that the approximate 
scattering matrix can then be written as a sum on poles under the form 



S (k) 



-2ika 



2N 



i + ikj2 



n=l k,Q ^ (fcg ^ /u) 

Equivalently, a product expression for the scattering matrix reads [Eq. @-(61)] 



(36) 



S (k) 



-2ika 



2N An) 



n 



w + k 



n=l rvQ Iv 



(37) 



This elegant result is valid only for the s wave. 

The internal wave function (r < a) is given by equation 
present notations 



2N An), 



u {r) 



-ike 



-2ika 



E 



'0 



1 />)fi» 

n=l n-Q V 



k) 



-(57) which reads in the 



(38) 



Equations ( |36D and (|38|) assume that the Siegert pseudostates (p^ are properly normalized 
[see Eq. @-(28)]. 



B. Equivalence for identical bases 



The above equations now allow us to prove the equivalence between the i?-matrix tech- 
nique of Refs. 0,0] and the Siegert-pseudostate method of Refs. . We shall first show 
that the approximations giving the scattering matrix are identical for any common finite 
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basis without mesh approximation. To this end, we specialize to I = and to short-range 
potentials. 

In Refs. PTT|JTm , the boundary parameter B is taken equal to zero but, as proved in 
Sec. Ill, exactly the same results would be obtained with any other value. Therefore we now 
focus on the choice B = Lq with the s-wave logarithmic derivative (|34]). The penetration 
factor Pq and the hard-sphere phase shift 0o take the simple forms 

P (ka) = ka (39) 

and 

4>o(ka) = ka. (40) 

Hence, Eq. (|T7p reads 

S (k) = e~ 2ika [1 + 2ikaR {ika)}. (41) 

This expression has the same structure as Eq. |fJ]-(58), since the R matrix is known to be 
related to the Green function through 

R l = {2a)- 1 G l {a,a), (42) 

see Eq. (IV. 1.10) of Ref. ||. It is thus equivalent to Eq. ( j37[) and relates the S'-matrix and 
complex i?-matrix poles. 

In order to prove the equivalence of both methods, we have to compare the approximate 
calculations of these expressions for finite bases. In Ref. [0], the Green function is obtained 
with Eq. ||-(49), where the matrix is obtained by inversion from Eq. ||-(50). Since the 
matrix appearing in Eq. j2|-(50) is identical to the present matrix C(L ) — EN [see Eq. (|8|)], 
the equivalence with our expression (|12"|) is proved. 



C. Equivalence of the mesh treatments 

We have just shown that, with the same finite basis, both methods are exactly equivalent. 
Now we show that the same property holds for the Discrete- Variable- Representation (DVR) 
approximation || of Ref. and the Lagrange-mesh approximation of Ref. \Y1\ which is 



summarized in section III. For / = 0, both methods make use of a mesh approximation 
related to zeros of Legendre polynomials. 



The R matrix is given by Eq. (|T|) with matrix C calculated with expression (|T9| ) where 



the different terms can be obtained on the mesh from Eqs. (|24|) and (|26| ) - (^3). In fact the 



DVR approximation for the Green-function matrix in Ref. |2j is the inverse of a matrix which 
is proportional to the Lagrange-mesh approximation of matrix Ciika) — EN [Eq. (|19~D1. 

The relations between the present quantities and those of Ref. are detailed in Table 
|. The first line of the Table shows that the notation x\ represents different zeros in both 
papers: in Ref. 0, they are zeros of a standard Legendre polynomial Pn(x) and belong 
to [—1, 1] while here and in Refs. |l~T| , p!2| they are zeros of a shifted Legendre polynomial 
P/v(2x — 1) and belong to [0, 1]. This is only a notational difference. The second line shows 
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the connection between the Gauss weights. The factor of two arises from the different lengths 
of the intervals. The basis functions are related in the third line. The functions TTi(x) in 
Ref. are chosen in such a way that they provide a representation of the unit operator 
[Eqs. |§-(5) and @-(6)] (see also Ref. [JD]). The property [0-(C13) of these functions shows 
that they verify a Lagrange condition. A compact expression for them has already been 
derived in Ref. || and used in Refs. |IT|Jl^| [see the present Eq. 

ha(x + 1) but, in Ref. 



1|)J. The present functions 
this factor is included in 



^ and used in Refs. 
r) may seem to differ by a factor r = ha(x + 1) but, in Ref. 
the operator as in Ref. [|ll|]. The overlaps differ from the present by a simple factor. 
Notice that the present form fl2"4| ) is simpler than expression 0-(C22) because we do not 
employ two different bases. 

Since the potential matrix elements are approximated in the same way, only the equiv- 
alence of the treatments of the kinetic energy remains to be proved. Let us detail the 
derivation of expressions ( p7| ) and (p8|) . Up to a factor 1/2, these matrix elements become 
with the Gauss quadrature 



(43) 



fi(r)fl(r)dr + J q /,(r)5(r - a)fi(r)dr 
= -{a\ i ) 1 ' 2 f'{ax i )+f i (a)fAa). 



The Gauss quadrature is exact for polynomials up to degree 2N — 1 |16[ so that expression 
( [13]) is exact and leads to Eqs. (|27| ) and fl2"8| ) (see Ref. |12| for technical details). Equivalently, 
Eq. ( |4~3"1) can be written as 



fi{r)f'{r)dr 



N 

a D ^kfl{axk)fj(ax k ). 

k=l 



(44) 



This expression is also exact but less compact. Still another approach is used in Ref. 0: the 
left-hand side of Eq. (H4T) is evaluated analytically after expressing the functions /j in the 
basis of Legendre polynomials, in the spirit of the DVR method |4]] . The notations for the 
kinetic-energy matrix elements are compared in the fifth line of Table Q. Comparing Eqs. ||- 
(C20) and ||-(C21) for Kij with the present Eqs. (^) and (|28| ) shows that the kinetic-energy 
matrix elements are calculated much more easily in the Lagrange-mesh philosophy than in 
the DVR philosophy. We have checked numerically that the expression |2|]-(C20) provides 
exactly the same results as ours, as it should. The notations for the full matrices are 
compared in the last line of Table |. 

Finally, let us mention a difference between the practical applications of the methods 

ll|,12|. In Refs. |iT,T2|, expression 



of Ref. M and of Refs. 



has been replaced by its 
the overlap matrix N is replaced by the unit matrix. This 
In Ref. ||, the algorithm requires that the equivalent 

This can easily be 



Gauss approximation (25): 
simplification is not used in Ref 

of matrix Ciika) be multiplied to the left and to the right by iV -1 / 2 
done in the present framework (see Sec.V) with iV~ 1//2 calculated as explained in Appendix 
A. However this complication is useless at the practical level because it does not improve 
the accuracy [0. We shall come back on the interest of the simplification (^) when dealing 
with the first example in Sec. V. 
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D. Consequences 



After those lengthy but necessary technical considerations, let us try to learn some 
practical consequences from the equivalence of the methods. 

For s states, the equivalence of both methods provides a new approach to the determi- 
nation of the poles of the S matrix in the i?-matrix formalism. Indeed deriving the complex 
S'-matrix poles from the real i?-matrix poles is not obvious. However the Siegert pseu- 
dostates (f)^ appear naturally in the .R-matrix formalism. They are solutions of Eq. fl3~T|) 
with complex wave numbers k\ n \ For I = 0, the algebraic algorithm of Ref. |2| provides an 
efficient way of determining some poles of the S'-matrix with sufficient accuracy. Note that 
only a few physical poles need usually be determined since the S matrix is more conveniently 
given by Eq. (0) than by Eqs. © or fl37|). 

Symmetrically the equivalence is also useful to attack the same problem for higher partial 
waves. The authors of Ref. have tried without success to generalize their search of 
Siegert pseudostates to / > 0. They make use of Jacobi polynomials adapted to the value 
of I. This basis was shown in Ref. [HJ to be not more efficient than basis ([H]), but more 



complicated to use. The reason for this failure is not due to a technical choice of basis 
but seems rather rooted in the technique of calculation of the scattering matrix. Indeed, as 
illustrated later in the second example below, the approximate wave functions determined 
in Ref. should be accurate, up to a normalization factor. As shown at the end of Sec. II, 
the same approximate wave functions are obtained for any choice of boundary parameter, 
independently of its physical adequacy. This includes the complex choice (|4]) which is 
implicit in Ref. 0. Hence only the scattering matrix is inaccurate in Ref. ||. 

The natural generalization of the Siegert-pseudostate method to I > and long-range 
potentials is Eq. (|33|) . However this equation does not allow to use the algebraic technique 
because the non linearity is not any more quadratic. We think that it would be useful to use 
Eq. ( |33| ) anyway to derive physical poles of the S matrix. The fact that it seems hopeless to 
find in this way all the pseudostates is of little importance since the S matrix can easily and 
accurately be calculated with the i?-matrix equation ([13]). The search for the physical poles 
could for example be performed by extending the iterative algorithm of Descouvemont and 
Vincke [ffHI. 



V. ILLUSTRATIVE EXAMPLES 

A. Bargmann potential 

Many examples are treated in Ref. @ and we have reproduced these results. In the 
case of phase shifts, we have checked that we obtain the same values both with the Siegert- 
pseudostate method and with the i?-matrix method within the accuracies of both numerical 
algorithms. When the number N of mesh points is not large enough, both methods provide 
essentially identical inaccurate results. Rather than repeating here one of those examples, 
we have chosen a different one which provides interesting intuitive information on the notion 
of Siegert pseudostate. We also use this example to discuss more deeply the effect of the 
Gauss approximation fl25|) on the overlap matrix element (]24f) . 
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The Bargmann potential |T9| , |2~0"H is denned as 



-2br 

V (r) = + (45) 

with (3 = (b — c)/(b + c) where b and c are real parameters. This potential has the remarkable 
property that its Jost function has only one pole and one zero 

- (46) 

The potential has one bound state for c < or one virtual state for c > 0. The scattering 
matrix reads 

So[k) ~ W ~ (k-m + icy (47) 

It possesses the symmetry property S (bc/k) = S (k). How approximation fl57p simulates 
expression (|i7D is instructive. 

The wave functions of the Siegert states of the potential truncated at a read 

where feg"' 1 is a solution of the equation 



b 2 (c + ik) . , 

5 sm ka = k 

cosh 6a 



(6 2 + zA;c) tanh 2 fca + b(c + i/c) tanh ba 
k 2 + b 2 



b 2 — c : 



(b + ctanh ba) z 



exp(— i/ca). (49) 



The unique Siegert state of the Bargmann potential fl45|) is given by Eq. ( ]48|) with the wave 
number — ic which is an approximate solution of (f49"D when a is large. 

The values 6 = 2 and c = — 1 are selected, which lead to a single bound state at energy 
— 1/2. Exact wave numbers of the Bargmann potential truncated at a = 5, i.e. replaced by 
zero beyond that value, are displayed as dots in Fig. [I] and in the first two columns of Table 
|I[ The isolated wave number close to i corresponds to the bound state. It is not exactly i 
because of the finite value of a. The pole of the simple Jost function ( |4*B| ) is simulated by 
a line of wave numbers with imaginary parts close to — 2i and rather regularly spaced real 
parts. The accuracies of the different approximations can be checked with respect to these 
exact values. 

In a first step, we perform a calculation of the wave numbers corresponding to the 
Siegert pseudostates by following the technique of Ref. without further approximations. 
For a = 5 and N = 25, the obtained values are depicted as circles in Fig. |l|. For real parts 
comprised between about —6.5 and 6.5, they closely correspond to the exact wave numbers 
of the cutoff Bargmann potential [see panel (a) of Fig. |l[] . The remaining poles split from the 
exact ones. They do not have any physical significance. They correspond to unconverged 
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solutions of Eq. fl35|) as shown by the comparison with the exact values in Fig. |l|. When N 
is increased, the horizontal line of physical poles of the cutoff potential obtained with the 
Siegert-pseudostate method extends to larger values and the non physical branches move 
accordingly. 

Numerical values are provided in Table [TT] for N = 25. The first ones agree within 
better than 10~ 5 with the exact values. These results are unstable: tiny changes in the 
computational algorithm (or a change of computer) modify the non converged digits. This 
instability appears to be related to factors 1 — X{ in the denominators of expressions (p7D . 
(P8Q, and fl3D|), for zeros Xj close to 1. We have also observed that the accuracy on the first 
wave numbers progressively deteriorates when N becomes larger than 25. This behavior is 
due to the fact that the norms of the corresponding eigenvectors are very small so that few 
digits of the eigenvalues are significant. When N increases beyond some value, the number of 
significant digits of each wave number decreases but the number of physically significant wave 
numbers nevertheless increases as mentioned above. Because of this important instability of 
the results for large N, the value N = 25 is close to optimal for a = 5 in the present case. 

Now we add a simplification to the method by replacing the overlap matrix A?" by a unit 
matrix according to the Gauss approximation fl25|) . This approximation has been used in the 
.R-matrix treatment of Refs. ]T2| , P^| , but it does not introduce much simplification in that 



case. The algorithm of Ref. [g] is much more simplified by this approximation. The obtained 
wave numbers are depicted as crosses in Fig. [I]. All the physical poles of the cutoff potential 
remain essentially unmodified (see details in Table [□]). On the contrary, the structure of the 
unphysical poles is completely different [see panel (b) of Fig. []J. As we shall see, this does 
not affect the phase shift at low energies. 

The first values obtained with approximation (25) are essentially identical to those with- 
out that approximation (see Table |TJ). This is very surprising because the agreement is much 
better than the accuracy on these values. When the real part increases, the values become 
more different, as already illustrated by the figure, but the results obtained with the Gauss 
approximation first remain better than those obtained with the exact overlap matrix [see 
panel (a) of Fig. [I]. This effect is also observed in other applications of the Lagrange-mesh 
method [|K)J. 

Phase shifts calculated for different energies with Eq. ( |37D are displayed in Table |T| and 
compared with the exact values from Eq. (|47|). Because of the symmetry property of the S 
matrix, the phase shifts at energies E and b 2 c 2 /4E add to ir. The calculations are performed 
for two choices of channel radius a, i.e. 5 and 6, and for two numbers of mesh points, i.e. 25 
and 40. 

For a = 5 and iV = 25, the relative accuracy on the phase shift is about 10~ 7 , which is 
better than for the wave numbers. We have also calculated the phase shifts with Eq. (p6|). 
The results are inaccurate as mentioned in Ref. ||] because the normalization of the Siegert 
pseudostates is difficult to achieve numerically: both terms of the normalization formula ||- 
(26) nearly cancel each other for several states. In the following we only refer to Eq. (j37[). 
Essentially the same values are obtained with the .R-matrix mesh method and with both 
variants of the Siegert-pseudostate method. A difference only appears at E = 10 but is not 
significant. The three results are much closer to each other than to the exact one. Increasing 
the number of mesh points to N = 40 without modifying a does not really improve the 
situation at E = 0.1 but is more useful at higher energies. The remaining disagreement is 
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due to the value of a for which the potential is not yet fully negligible. 

For a = 6 and N = 25, the relative accuracy on the phase shift is only about 10~ 5 . The 
number of mesh points is too small for this a value. At all energies, the three approximate 
results are much closer to each other than to the exact phase shift. When N is increased to 
40, the relative accuracy on the phase shifts is better than 10 -9 . No significant differences 
between the methods appear. The influence of the approximation (p5|) is weak for all the 
considered energies. 



B. Purely centrifugal potential 

In Ref. |§, the Siegert-pseudostate method is applied to a purely centrifugal potential. 
In this case, the natural choice ( |i~6|) for B would be 

since 0\(kr) = — 1/ikr) exp(ikr). The choice B = L = ika is made instead in order 
to use the algebraic algorithm established for I = short-range potentials. The S matrix is 
then calculated with Eq. (|^) which leads to very inaccurate results as can be expected in 
the present framework. 

The corresponding wave function is however essentially independent of the particular 
choice for B (see Sec. II). Hence, it is as accurate for any non-physical value of B as for 
the physical one ([5(]). However, its normalization is not correct. We have verified that the 
method of Ref. @ and in particular Eq. [@]-(57) or the present Eq. fl38|) provide accurate 
free-particle wave functions for the centrifugal potential, up to a multiplicative factor. This 
fact seems to have remained unnoticed by the authors of Ref. . It can be understood from 
the present Eq. ( [T5| ) to which Eq. Q58| ) is equivalent. Expression (|T5|) is valid for any B but 
the correct normalization factor Ii(ka) — SiOiika) is not available to the Siegert-pseudostate 
method since L is used instead of Li and the calculated S matrix is inaccurate. Moreover, 
this correct scattering wave function is obtained with unphysical Siegert pseudostates and 
S matrix poles. These unphysical states are just used as a basis for the expansion of the 
wave function. 

Our result has been obtained with basis (|2lf) based on Legendre polynomials without 
recourse to the more complicated Jacobi polynomials employed in Ref. 0. This confirms 
the fact established in Ref. |12| that the /-dependent basis of Jacobi polynomials, although 
accurate, is not necessary to treat I > partial waves. 



VI. CONCLUSIONS 

In this paper, we have proved that the Siegert-pseudostate and Lagrange-mesh i?-matrix 
methods are in fact exactly equivalent. This property is true when the same finite basis 
is used in both approaches but also when the mesh methods described in Ref. || and in 
Refs. (12|;[13| are employed. 

This equivalence provides an approximate way of calculating the poles of the S matrix in 
the i?-matrix framework for s wave short-range potentials. It also shows how to generalize 
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the determination of Siegert states for / > and for long-range potentials but the algebraic 
algorithm developed in Ref. || must be replaced, most probably by some iterative technique. 

We have shown that, unexpectedly, the method of Ref. [0 can also be used to construct 
scattering wave functions for potentials with long-range terms such as the centrifugal term. 
However, in such cases, the i?-matrix method of Refs. [|12| , |i"3f is more advantageous since 
(i) it is simpler, (ii) it also provides the S matrix, and (iii) it is readily extended to the 
multichannel case. 
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APPENDIX A 

Let B be an invertible matrix and u and v be vectors. The inverse of the matrix 

A = B + uv T (Al) 

is given by 

b-Vb- ^ 

1 + v T B' 1 u 

where the denominator is a scalar. A corollary of Eq. ( |A2| ) reads 

B~ x u 
1 +v T B~ L u 
Another corollary is the relation 

(v T A- l u)- 1 = 1 + (v T B- l u)- 1 (A4) 

from which Eq. (|1~4]) follows. 

The norm matrix has the form 

N = 1 + auu T . (A5) 

where u is here a unit vector (||m|| = 1). Arbitrary powers of N are given by 

N x = 1 + [(1 + a) x - l}uu T , (A6) 



A ~' u = ^r^rr- ( A3 ) 



for any A, integer or fractional, positive or negative. In Eq. (|24]), the components of the unit 
vector u read 



- - ^if-^r 1 (A7) 

and the coefficent a is given by 

a = N 2 /(2N +1). (A8) 
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FIG. 1. Complex wave numbers of the Siegert pseudostates of the Bargmann potential: exact 
(dots), without (circles) and with (crosses) Gauss approximation for the overlap matrix (a = 5 and 
AT = 25). 
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TABLES 

TABLE I. Symbols of Ref. @ expressed in the present notations (1 = 0). 



Ref. § 

X { 
Ki 
TTi(x) 
Pij 

K, 
L 



H + (1 - iJfco)X 



Eq. 

~W~ 

(Cll) 
(C9) 
(C19) 
(C20) 
(C18) 
(C15) 



Present 
2xi-l 

to)V2 Xi ( x + i)-i/.[i a ( x + i)] 
a %i%j{fi\fj) — XiXjNij 

\a 2 x i Xj{f i \r~ 1 5(r - a)\fj) 
a 2 XiXjC(ika) 



Eq. 

w 
m 

(H 

(1) 

(HI) 



TABLE II. Complex wave numbers of the Siegert pseudostates of the Bargmann potential: 
exact [Eq. (f49|)1, without [Eq. ([2^)1 and with [Eq. (p5|)] Gauss approximation for the overlap matrix 
(a = 5 and iV = 25). 





Exact 


Eq. ( 


H 


Eq. (25) 


Re k^ 


Im k^ 


Re k^ 


Im /cq^ 


Re 4 n) 


Im 4 n) 





0.9999999999955 


0.0000000 


0.9999999 


0.0000000 


0.9999999 


0.6390266702 


-1.9909403640 


0.6390286 


-1.9909403 


0.6390286 


-1.9909403 


1.2947020189 


-1.9808654900 


1.2947061 


-1.9808653 


1.2947061 


-1.9808653 


1.9559919327 


-1.9824365206 


1.9559981 


-1.9824363 


1.9559981 


-1.9824363 


2.6142071140 


-1.9926145823 


2.6142153 


-1.9926144 


2.6142153 


-1.9926144 


3.2675033086 


-2.0070370919 


3.2675133 


-2.0070370 


3.2675135 


-2.0070369 


3.9163165163 


-2.0231436137 


3.9163283 


-2.0231433 


3.9163285 


-2.0231433 


4.5614908509 


-2.0396289360 


4.5615045 


-2.0396283 


4.5615048 


-2.0396285 


5.2037975084 


-2.0558538335 


5.2038193 


-2.0558336 


5.2038125 


-2.0558513 


5.8438509338 


-2.0715187039 


5.8453967 


-2.0710367 


5.8438594 


-2.0715509 


6.4821222431 


-2.0864975209 


6.4948355 


-2.1966169 


6.4817804 


-2.0855307 


7.1189698354 


-2.1007533271 


6.6971747 


-2.0045040 


7.1321815 


-2.0820124 


7.7546673954 


-2.1142944575 


7.2971683 


-1.5698719 


7.5590343 


-2.4165539 
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TABLE III. Phase shifts (in degrees) of the Bargmann potential: exact [Eq. (f47|)], Siegert 
pseudostate method without [Eq. (|24|)] and with [Eq. (^)] Gauss approximation for the overlap 
matrix, and R matrix [Eq. fll3|)1. 



E 


Exact 


a 


N 


Eq. © 


Eq. (H) 


Eq. (ID 


0.1 


143.30077480 


5 


25 


143.30076489 


143.30076489 


143.30076490 








40 


143.30076866 


143.30076866 


143.30076866 






6 


25 


143.30079803 


143.30079803 


143.30079803 








40 


143.30077472 


143.30077472 


143.30077473 


1 


90.00000000 


5 


25 


89.99999083 


89.99999083 


89.99999083 








40 


89.99999947 


89.99999947 


89.99999958 






6 


25 


90.00005685 


90.00005678 


90.00005688 








40 


89.99999998 


89.99999998 


90.00000007 


10 


36.69922520 


5 


25 


36.69916732 


36.69916703 


36.69916730 








40 


36.69922475 


36.69922475 


36.69922502 






6 


25 


36.69975390 


36.69974764 


36.69974770 








40 


36.69922520 


36.69922520 


36.69922526 
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